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Abstract. A recent approach for constructing an exact 3D reduction of the 4D matrix 
elements of the electromagnetic current for an interacting two-fermion system, is briefly 
reviewed. The properties of the obtained 3D current, like the fulfillment of the Ward-Takashi 
Identity and its Fock decomposition, will be illustrated in relation with future applications to 
^ ' few-nucleon systems. 
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CN ' 1. Introduction 

A detailed investigation of the electromagnetic (em) properties of hadrons and nuclei requests a 
careful treatment of the current operator for interacting fermion systems. In particular, besides 
general properties like the Poincare covariance, the gauge symmetry of the em interaction plays 
the well-known essential role, that leads to the Ward-Takahashi Identity (WTI) for the em vertex 
function. In this short presentation, some insights will be given on the main ingredients of our 
approach [1] aimed to construct a 3D em current for an interacting two-fermion system, that i) 
fulfills WTI and therefore allows the current conservation (CC) and ii) preserves at large extent 
the Poincare covariance, emphasizing the relevant steps for obtaining workable expressions to 
be used in actual calculations. 

Field theoretical approaches, based on the Bethe-Salpeter equation, (see, e.g. [21 S] ) allow 
one to evaluate matrix elements of the 4D em current for an interacting two-body system, viz 

(P f ,J f ,a f \f\P i ,J u a i )=CY / J j^^yAhX^Pf^ih^^i^^Pi) (1) 

where, in the present case, ^ Jcr (k\, &2, P) is the Bethe-Salpeter (BS) amplitude of a two-fermion 
systems (with k± + k' 2 = Pt and k\ + k 2 = P% for the final and initial state, respectively), C 
contains the normalization factors, and is the 4D fermion-photon vertex (or kernel) that 
contains one + two +...+ many-body contributions, and must fulfill WTI. The sum in Eq. ([1]) 
is performed on isospin, spinorial, etc. indexes. 



To solve the BS equation in Minkowski space for obtaining (k\,k2, P), in realistic cases, 
is still a prohibitive task, but new and very effective approaches are coming into play (see, e.g. 
[5]). On the other hand, 3D reductions are well known tools (see, e.g. [6]) for attempting to 
overcome the above mentioned difficulties, but still saving at large extent the general properties 
that the em current has to fulfill. 

Aim of our investigation pQ is i) to achieve an exact 3D reduction of the 4D matrix elements 
in Eq. (pTJ) , i.e. without loosing any physics content, and ii) to elaborate a systematic analysis 
of the obtained expression in order to develop effective approximation. The main ingredients 
for our approach, that were already applied to the case of two-boson interacting system in 
[3], are in order i) the Quasi-Potential (QP) approach for the Transition matrix (see [7]), and 
ii) the projection of the 4D physical quantities onto the 3D Light-front (LF) hyperplane (i.e. 
x+ = x° + x 3 = 0) (see, e.g. [2J As a first result, one can establish a formally exact 
correspondence between 4D BS amplitude and the 3D LF "valence" wave function. Then it is 
possible to express, without approximation, the matrix elements of the 4D em current, where 
the BS amplitudes appear, in terms of matrix elements of a 3D LF current, that fulfills the WTI. 
It should be pointed out that in these last matrix elements the valence states are present. 

Within the Quasi-Potential approach for the T-matrix, where the so-called effective 
interaction appears, one can develop a meaningful approximation scheme for evaluating the 
3D matrix elements, based on the Fock decomposition of the effective interaction. In particular, 
one expects a decreasing influence of the Fock states that contain an increasing number of 
intermediate particles. This machinery allows one to construct an iterated approximation of 
the effective interaction and eventually to evaluate matrix elements of the "truncated" LF em 
current, that still fulfills at any order the LF WTI. Finally, by using a Yukawa model in ladder 
approximation, as a pedagogical example, explicit expressions for many-body contributions to 
the em LF current have been obtained pQ. 

2. BS amplitudes and EM current 

The quantities present in Eq. ([1]) are the BS amplitude and the em vertex. In particular, the 
BS amplitude \^/), that depends upon internal variables, satisfies the following equation (the 
lim e ^o is understood) 

G- 1 (K)\^) = (2) 

with appropriate boundary conditions for bound and scattering states. In Eq. ([2]), K is the 
4-momentum of the center of mass and the inverse of the Green function is given by 

G- 1 (K) = G 1 (K)-V(K) (3) 

where V(K) is the 4D interaction, as obtained by the Lagrangian, Gq(K) = S{k\) S(K — k\), 
with S(k) the Dirac propagator. 

The em vertex, J 11 , corresponding to the interaction V(K) has a free term, Jq , and another 
one, J i ', which depends on the interaction, as requested by the commutation rules with the 
generators of the Poincare group. It also satisfies the following Ward-Takahashi identity 

QtiJ^iQ) = G-\K f )e - e G" 1 ^) (4) 

where is the 4-momentum transfer and e = e± + ei is the pointlike charge operator with 
matrix elements given by (kj\ej\pj) = ej5 4 (kj — pj — Q) 

Combining Eqs. (J3j> and dH) with the decomposition of in terms of the dependence upon 
V(K), it turns out that each contribution fulfills a suitable constraint, viz 

QMQ) = G Q \K f )e - e G, 1 ^) 



and 

Q^Ji\Q) = eV{K i )-V{K f )e . 
Finally, CC can be explicitly retrieved by using solutions of the BS equation, Eq. ([2]), i.e., 

Qr{* J \J»{Q)\* i ) = (V f \ [G~\K f )e -eG- l (Ki)] = . (5) 

After introducing the general equations fulfilled by BS amplitudes and em current, let us 
illustrate our procedure. As well known, the T-matrix is given in terms of the interaction V(K ) 
by an integral equation: T{K) = V{K) + V(K) G (K) T(K). 

In the Quasi- Potential framework [7j, T{K) can be rewritten in terms of two coupled equations 

T(K) = W(K) + W(K) G aux (K) T(K) (6) 

where the effective interaction W(K), in turn, is a solution of 

W(K) = V{K) + V(K) A (K) W{K) (7) 

with Aq(K) := Gq{K) — G aux {K) , and G aux {K) is a smart choice to speed up the convergence 
of the iterative solutions. 

3. Projecting onto the Light-Front hyperplane 

In Eq. ([I]) the 4D integration contains all the difficulty given by the very complicated analytic 
structure of both BS amplitudes and em vertex in Minkowski space. The adopted strategy for 
decreasing the degree of difficulty, is based on disentangling the "dynamical" variable, i.e. the 
time, from the remaining three kinematical variables, that describe the "initial" state. In order 
to accomplish such a task, the LF kinematics appears a very appealing tool (see, e.g. [8] for 
an extensive review) since the LF combination of the Poincare generators has the property to 
produce the maximal numbers of generators without interaction (7 out of 10). In turn, such a 
property produces many other relevant features, like the kinematical nature of the LF boosts 
(to be not confused with the standard boosts, a part the longitudinal one). In particular, the 
LF components of the 4-momentum and the scalar product are given by 

k ± =k ±k Z , k ± = {k X ,ky}, k-X= ( X k+ + X+k l _ X± ■ k± . ( 8 ) 

Integrating a given physical quantity, (fi(k), over the minus component of the 4-momentum, k~ , 
amounts to restrict its dependence upon only 3 variables in the coordinate space. One has 

J dk~ <p(k) = J dk~ 1 J dx + dx~ dxj_ e lk ' x <p(x) = 



1 



(2tt) 2 



2(2tt) 4 

dx+ dx- 5(x + )d* ± e i(*-k+/2-x ± -k ± ) ^ _ (g) 



Therefore, the "LF-time", x + , that labels the dynamical evolution (as the standard "t" in the 
instant form of the relativistic Hamiltonian) , is constrained to its initial value, i.e. x + = 0. 
Therefore, once the physical quantities are integrated over k~ , they are constrained to live onto 
the LF hyperplane defined by x + = 0. 

A crucial step for the LF projection is the separation of the fermion propagator in an on-shell 
term and in an instantaneous one, i.e 

% S{k) = T — = — — — j— + — — = i S on (k) + — — (10) 

k 2 -m 2 + ie h+ (k~ - k on + ik) 2fc+ 2k+ 



where k~ n = (k^_ + m 2 )/k + is the minus-component of k£ n , such that k on • k on = m 2 , and the 
second term leads to an instantaneous (in LF time!) free propagation, since the formal Fourier 
transform generates 8{x + ) (see Eq. ©). It should be pointed out that such a term makes the 
treatment of a fermionic system basically different from the treatment of a bosonic one. On the 
other hand, the decomposition in Eq. ()10p suggests a strategy for achieving our goal of an exact 
3D reduction. As a matter of fact, let us define the 4D on-shell two-body Green's function [HH] 

G° n (K) := S on (h) S on (K - h) (11) 

that plays a key role in our approach, since allows us to separate the trivial 4D propagation 
associated to such a contribution. 

The 3D LF counterpart of Go™ (10 is introduced by integrating over k^ and k'{~ , i.e. by 
projecting onto the 3D LF hyperplane, the following matrix elements (note that the other three 
components of the momentum remain operators) 

go(lO = lG° n (K)l = n+ j dk'fdki {k' x ~\G^{K)\k^) (12) 

where Q + is a suitable kinematical factor ensuring the kinematical Poincare covariance (e.g. 
Q+ = k + or Q + = ^Jk + {K + — k + )), and the symbol 1, from now on, indicates the integration 
over the minus component of the 4-momentum that appears on the left or on the right side, 
depending upon the position of the symbol itself. 

In the case without interaction, it is easy to find the relation between the 3D LF free state, 
|0o), solution of 0(5~ 1 (1C) |0q) = (note that instantaneous terms do not contribute) and the free 
BS amplitude, solution of Gq 1 (K) |\&o) = (note that the particles are on their own mass-shell 
in the free case). It turns out pQ that 

|tt„> = G% n (K)l g^\K) |0 O ) = IL (K) |</> ) (13) 

where Ho (10 is the free reverse LF projector, that allows to jump from a 3D description to a 
4D world. The operator Ho (10 has a simple explicit expression [TJ E] in terms of the positive 
energy Dirac projector. Let us remind that, by construction ^l^o) = |^o) : this completes the 
reconstruction pattern in the simple case of a free two-fermion system. 

In order to apply the projection technique to the QP approach, namely for an interacting 
case, one needs an intermediate step, given by a suitable choice for the auxiliary Green's function 
to be used in the 4D operator Ao, see Eq. ([7]). The chosen 4D auxiliary Green's function is 
defined as follows in terms of the 3D go{K) 

G aux {k) = U (K)g (K)U (K) (14) 

where Ho (10 = ^^{K). Then, the T- matrix and the effective interaction W(K) can 

be rewritten as follows 

t(k) = w(k) + w(k) [n (i0<7o(i0no(i0] T{K) 

W{K) = V{K) + V(K) A (K) W{K) = 

= V(K) + V(K) [g (K) - (n (K)g (K)U (K))] W(K) . (15) 
The equation for W(K) can be solved by iteration, i.e. 



oo 

w(k) = v{k) j2 [Aownior 1 



(16) 



The rate of convergence of this series is related to Aq(K), that describes the difference between 
the intermediate propagation of two non-interacting fermions and the 4D counterpart of the LF 
propagator g${K). Depending upon the powers of Aq(K), we will have intermediate states with 
more and more bosons and fermions (once we introduce time-ordered diagrams). 

The Rosetta Stone relating 4D and 3D quantities, is readily constructed. The 3D T-matrix 
and effective interaction are respectively 

t(K) = U (K) T(K) U (K) w(K) = U {K) W{K) U (K) . (17) 

Once t(K) is reexpressed in terms of go{K) and w(K), one can introduce the 3D interacting 
Green's function and its inverse as follows 

g(K) = g (K) + g (K)t(K)g (K) g~\K) = U Q (K) G~\K) U(K) (18) 



where the interacting LF reverse projection operator is 

U(K) = [1 + A (K) W{K) ) U (K) . (19) 
The solution \<p) of the following 3D equation, 

g~\K) \<f>) = [g \K) - w{K)] \4>) = (20) 

with appropriate boundary conditions for bound and scattering states, is the valence component 
of the LF wave function for two interacting fermions, as can be deduced by the relation with the 
BS amplitude, shown in what follows. The exact relations between the 3D valence component 
and the 4D BS amplitude can be readily obtained, and they are given by 

m = IL(K) \4>) lG° n {K) Gq\K) |*> = |0) . (21) 

The rightmost equation illustrates the motivation for calling \<p) valence component, given the 
presence of the integration on the minus component in lhs (see, e.g [5] and references quoted 
therein). Such an integration makes only the first component of the Fock decomposition of \^>) 
to survive. 

A possible approximation scheme is suggested by the QP expansion of the effective interaction, 
Eq. (|16p . Then, w(K), t{K) and the 3D LF Green's function g(K) (all containing W{K) in 
a more or less explicit way) can be evaluated by truncating the expansion at a given order. It 
should be pointed out that the full complexity of the Fock-space affects g{K) through the effective 
interaction w(K), as can be seen through a time-ordered analysis of the diagrams associated to 
each contribution [Aq(K)V \K)\ l ~ , and reminding that Aq(K) describes propagation. Indeed, 
the truncation of the QP expansion limits the number of Fock states involved in the construction 
of the effective interaction. As a final remark, one could argue that the convergence rate of the 
QP expansion is related to the small probability of the higher Fock-components. 



4. The 3D EM current and the LF Ward-Takahashi Identity 

Once we have obtained the relation between the 4D BS amplitude and the 3D LF valence wave 
function, one can establish a direct link between the matrix element of the 4D current and the 
matrix elements of the 3D LF em current operator. 

For both scattering and bound states one has [U [U El [3] 



^ } \J^{Q)\m i ) = (4> f \f{Kf,K i )\ ( t> i ) 



(22) 



where the 3D LF current operator, acting on the valence wave functions, is defined as follows 



f(K f ,Ki) = U(K f )J^Q)U(Ki) = 

= U (K f ) [1 + W(K f )A (K f )] J»{Q) [1 + Ao(Ki)W(Ki)] n (i^ 



(23) 



Since the 4D current must fulfill the WTI, Eq. (jj]), the 3D LF current fulfills the following 
LF WTI PQ 



Q^f(K } ,K % ) = U(K f ) [G-^K^e-eG-^Ki)] 11(1$ 
= 9-\K } ) Ql F -Ql F g-\ Ki ) 



(24) 



1 ( R) 

where Q LF is called the left (right) LF charge operator (note that it is interaction free !). For 
the sake of concreteness, let us show the explicit expression of the left charge, viz 



mi 



QlLF = A+(fclon)-f— 7l ^ lLF K + (ki on ) A H 
Ki 



(25) 



where A + is the positive energy Dirac projector, and e.\LF the 3D pointlike charge operator, 
given by 



(k{',k 1 _ L \e 1LF \k{,k lj± ) := e x 5 lk'{ - k{ 



5* (k' 1± 



Q. 



(26) 



The current conservation straightforwardly follows by taking the matrix elements of Eq. (|27 
between 3D interacting states solutions of g~ l (K) \<f>) = 

For the actual calculations, one needs to truncate the iterated expression of the effective 
interaction W(K), Eq. (fl~6|) . and therefore one ends up with a truncated 3D em current. The 
truncation scheme for the LF current cannot be carried out in a naive way, if one has to fulfill 
a WTI. Indeed, at each truncation order it will correspond a truncated effective interaction, 
w( n \K), a truncated LF Green's function g( n \K), a truncated valence wave function, a 
properly truncated 3D LF current, and finally a truncated WTI. Let us see in details this chain. 

Once we truncate the QP expansion of the effective interaction at the order n, 



W (n \K) = ViK^lMKWiK)}*- 1 



(27) 



i=i 



we can immediately generate an approximate 3D LF Green's function, g( n \K), leading to an 
approximate valence wave functions, \(f>^ n '), solutions of the following equation 



[g^(K)]- 1 = \g,\K) - w^(K)} \<j>^) = 



(28) 



The corresponding em current, that fulfills CC, if one uses \(f)^) for evaluating the matrix 
elements, cannot be directly constructed inserting in the definition of the 3D current, Eq. 
the truncated W^ n \K), i.e. 



Q ■ f n \K f ,Ki) + [g^iKf)}- 1 Q L LF - Q* F \g^\Ki 



(29) 



with 



(K f , Ki) = U (K f ) \l + W^\K f )A (K f )] J»(Q) 



l + AoOfQ)W(")(iQ)l UoiKi 



(30) 



The fulfillment of the LF WTI excludes such a naive approximation scheme for the truncated 
current, but, at the same time, gives us the correct hint. One should implement a power counting 
of the effective interaction that appears in the definition of the LF current (the interaction V(K ) 
is present both in the projectors H(K) and in the 4D current J) and in the rhs of the LF WTI, 
where only [g^^Kf)] -1 and [</ n )(-fQ)] _1 are affected by the interaction, as in the exact case, 
i.e. Eq. USD. 

For the truncation at the n-th order, a truncated current that satisfies a truncated LF WTI 
is given by [T] 



c(n)/i _ jc(n-l)fi _|_ jj c 



n n— 1 

i=0 i=0 



n (3i) 



where Wi = V [AoVf 1 and it has been formally defined WoAo = AqWo = 1- It should be 
reminded that J$ is O(V ) and Jf is 0{V l ). 
Then j c ( n )^ satisfies a LF WTI given by 

Q% (n) = [g {n) }-HKf) q l lf - Qlf [g^r'm . (32) 

For obtaining CC, the matrix elements should be taken between solutions of [g ijl ^]~ 1 \<pn) = 0. 

5. Two-body current in the LF approach: a pedagogical example 

A Yukawa model in ladder approximation has been adopted in order to elaborate an application 
of our general approach. Within such a model, a couple of fermions interacts by the exchange 
of a boson, in the present case chargeless, for simplicity. From the Lagrangian, one can deduce 
immediately the 4D interaction to be used in the definition of the T-matrix and of the 4D 
current. In ladder approximation, the 4D current is nothing else but the free 4D current that 
fulfills the suitable WTI (see, e.g. [9]). Let us stress that, though in Minkowski space the current 
of the model under investigation is the free one, its projection onto a 3D hyperplane becomes 
interaction dependent. The matrix elements of the first-order LF current, in a 3D Fourier space, 
are given by 

k[ ± \f^\kfk 1± ) = {k'+k' 1± \ [f + TTo [V A0J5" + Jo"A V] Uo] \kfk 1± ) (33) 



','+2 



with j c (°^ = n J^U Q (remind that n is the free -case projector). In figure 1, the diagrammatic 
analysis of the first-order LF current is shown. It should be pointed out that i) the diagrams 
are LF-time ordered, ii) the LF-time flows from the right to the left and iii) the external 
legs are particles on their-own mass shell (therefore only three components can independently 
change). Such a diagrammatic analysis allows one to appreciate the physical content inside 

c(l) 

the operator . The peculiar feature of the LF approach is the presence of instantaneous 



contributions, generated by the term proportional to 7+ in Eq. (|10j) . This contributions 
are labeled by horizontal dashes. Another interesting contribution is the Z-diagram, or pair 
production diagram, generated by the time-ordering of the 4D Dirac propagator. In this respect, 
we should emphasize that the full Dirac structure, i.e. the spinorial structure, is exactly taken 
into account, since the LF projection only works on the minus components of the four-momenta. 

This first-order current will be applied for evaluating the Deuteron electromagnetic properties, 
in order to extend the analysis at the zero-order already performed [10]. But, one can anticipate 
that i) the pair term affects all the three electromagnetic form factors, while the instantaneous 
terms contribute only to the magnetic one (since the property 7 + 7 + = is acting in the 
charge form factor), ii) the pair term vanishes for q + — > 0, for the conservation law of the 
plus components, while the instantaneous ones survive, iii) the pair term should be maximal at 
q + ~ in 7v [11], iv) the remaining terms affect all the Deuteron form factors, in the whole range 




Kf — k[ K t — kx Kf — k[ K, — fc'i Kf — k[ K t — ki 

Figure 1. LF-time ordered diagrams for the first-order LF current for a Yukawa model in 
the ladder approximation, see Eq. ([33]) and text. Horizontal dashes indicate instantaneous 
propagations. (Adapted from [1]) 

6. Conclusions 

For two interacting fermions, an exact correspondence between 4D BS amplitudes and 3D LF 
valence wave functions has been obtained [1], by using the projection onto the LF hyperplane. 
The approach fulfills the covariance with respect to the kinematical Poincare subgroup (7 
generators out of 10). 

The obtained relation allows one to express the matrix elements of the 4D em current in 
terms of matrix elements of the 3D LF em current between valence wave functions. For such a 
current operator, the corresponding LF WTI can be constructed by introducing interaction free 
LF charge operators, then current conservation can be trivially retrieved. 

Finally, by using a suitable auxiliary Green's function, within the Quasi-Potential approach, 
that allows an ordering in terms of intermediate Fock states, an approximation scheme has been 
developed, for both dynamical equation and the em current operator. In particular, at any order 
of the effective interaction the em current operator fulfills a proper LF WTI, leading to CC. 

A systematic analysis of LF two-body currents, obtained within a Yukawa model in ladder 
approximation, is in progress for the Deuteron case. 
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